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We present quantitative measurements and calculations of the spin-orbit induced zero-magnetic- 
field spin-splitting in two-dimensional (2D) hole systems in modulation-doped GaAs (3f f )A quantum 
wells. The results show that the splitting is large and tunable. In particular, via a combination of 
back- and front-gate biases, we can tune the splitting while keeping the 2D hole density constant. 
The data also reveal a surprising result regarding the magnetoresistance (Shubnikov-de Haas) oscil- 
lations in a 2D system with spin-split energy bands: the frequencies of the oscillations are not simply 
related to the population of the spin-subbands. Next we concentrate on the metallic-like behavior 
observed in these 2D holes and its relation to spin-splitting. The data indicate that the metallic 
behavior is more pronounced when two spin-subbands with unequal populations are occupied. Our 
measurements of the magnetoresistance of these 2D hole systems with an in-plane magnetic field 
corroborate this conclusion: while the system is metallic at zero magnetic field, it turns insulating 
when one of the spin-subbands is depopulated at high magnetic field. 
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I. INTRODUCTION 



Zero-magnetic-field (B = 0) spin-splitting of energy 
bands has long been a subject of considerable experimen- 
tal and theoretical effort because of its role in the elec- 
tronic properties of two-dimensional (2D) systems |lp|§. 
It concerns details of energy band structure that are of 
fundamental interest. The eigenstates in a system that 
has inversion symmetry in both space and time are spin- 
degenerate. If the system is made spatially inversion 
asymmetric, this degeneracy is lifted by spin-orbit in- 
teraction for all momentum k ^ 0, even in the absence of 
an external magnetic field. In 2D systems, we can create 
this spatial inversion asymmetry simply by making the 
potential which confines the carriers to 2D asymmetric. 

In this work, we study the B = spin-splitting in GaAs 
2D holes confined to a quantum well (QW) grown on the 
(311)A surface of a GaAs substrate. Section || provides 
some experimental details. In Section II], we show that 
by using two gates on the back and front of the sam- 
ple to control the electric field applied perpendicular to 
the QW, the asymmetry and therefore the spin-splitting 
can be tuned while the carrier density is kept constant. 
In the process we discover that the low-i? Shubnikov-de 
Haas oscillations are not as simply related to the B = 
spin-splitting as has often been presumed Q. 

Tuning the spin-splitting without changing the density 
allows us to isolate the effect of spin-splitting on various 
phenomena. In recent years there has been significant in- 
terest in the surprising metallic behavior of various high- 



In Section IV we show the 



quality 2D systems 
effect of spin-splitting on the metallic behavior of the 2D 
holes in GaAs @,||. 

Finally, motivated by recent theoretical predictions 
pl| , p2| , we have applied an in-plane B to these sam- 
ples and studied the magnetoresistance (Section |y|). We 
have found that, remarkably, the in-plane B depopulates 
the upper spin-subband at a rate which depends on the 
relative orientations of B and the crystal axes! This is 
a reflection of the anisotropic band structure and Zee- 
man splitting in the system. Furthermore, when the up- 
per spin-subband is depopulated, the metallic behavior 
turns to insulating behavior. We also see an in-plane 
magnetoresistance anisotropy that depends on the rela- 
tive orientations of B and the current /; the finite layer 
thickness of our 2D system may be responsible for this 
observed anisotropy p2|. 



II. EXPERIMENTAL 

A schematic cross section of our samples is shown in 
Fig. ^A. The samples are grown, via molecular beam 
epitaxy (MBE), on GaAs (311)A substrates. Each sam- 
ple contains a 20 nm-wide GaAs QW flanked by AlGaAs 
barrier layers which are modulation doped with Si. On 
GaAs (311)A substrates, Si is normally incorporated as 
an acceptor, leading to a high-quality 2D hole system in 
the GaAs QW. We make ohmic contacts to the 2D holes 
using In:Zn and alloying in a reducing atmosphere. As 
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shown in Fig. [l]A, our samples have back and front gates 
to change the density of the 2D hole system and also 
tune the symmetry of the confining QW potential. The 
back gate is made by contacting a layer of In left over on 
the back of the substrate from the MBE growth, and the 
front gate is evaporated metal (Ti-Au). 

GaAs 2D systems grown on (311)A substrates have an 
intrinsic mobility anisotropy due to an anisotropic inter- 
face roughness |2^j24|]. The mobility for current along 
the [233] direction can be up to several times the mobil- 
ity for current along the [Oil] direction. For this reason, 
we did our measurements on L-shaped Hall bars, which 
allow simultaneous measurement of the resistivity p for 
current along these two directions (Fig. |I]b). 

Measurements were done in dilution and 3 He refriger- 
ators, with base temperatures of 30 mK and 0.3 K re- 
spectively. The systems are fitted in a superconducting 
magnet capable of magnetic fields up to 16 T. The 3 He 
refrigerator has a tilting stage for the sample that can be 
rotated in-situ, so the relative angle between the mag- 
netic field and 2D hole system plane can be varied. 

Figure |^ shows an example of data for one of the sam- 
ples from both mobility directions, measured simultane- 
ously at 30 mK . The inset to A shows the \ow-B data. 
The main part of A demonstrates the high quality of the 
sample. The inset to B shows raw T-dependence of p 
data, and the main part of B shows it scaled to po, the 
resistivity at 30 mK. 



III. TUNING SPIN-SPLITTING 

We use the gates to tune the asymmetry of the QW by 
applying an electric field (E±) perpendicular to its plane. 
Figure [l]C schematically demonstrates the procedure. We 
set the front gate (Vfg) and back gate (Vbg) voltages, 
and measure p as a function of B. Then, at a small B, 
Vfg is increased and the change in the hole density is 
measured from the change in the Hall coefficient. Vbg 
is then reduced to recover the original density. This pro- 
cedure changes Ej_ while maintaining the same density 
to within 3%, and allows calculation of the change in E± 
from the way the gates affect the density. These steps are 
repeated until we have probed the range of Vfg an d Vbg 
that are accessible without causing gate leakage jl2| . In- 
creasing the magnitude of E± increases the asymmetry 
of the sample, which increases the spin-splitting. 

Figure ||A shows examples of the Shubnikov-de Haas 
(SdH) oscillations measured for current along the [Oil] 
direction in our 20 nm QW, following the procedure de- 
scribed above to change E± from about 5 kV/cm (point- 
ing towards the front gate) in the top panel through to 
about -6 kV/cm in the bottom panel. Beating can be 
clearly seen in the oscillations of all traces except the 
center trace where E± ~ 0. Fast Fourier transforms 
(FFTs) of p vs. B^ 1 quantify the frequencies present. 



Figure ||B shows the FFTs of the SdH oscillations at all 
of the measured sets of Vfg and Vbg- These frequencies 
have long been thought to be directly proportional to the 
spin-subband densities (p±) following 



e f ± 

P± = frJSdH- 



(1) 



With this interpretation it is clear in Fig. gB that we 
are tuning the spin-splitting through a minimum as we 
change E± from 5 kV/cm to -6 kV/cm: the two peaks, 
corresponding to two spin-subband densities, get closer 
together as E±_ approaches zero, finally merging, only to 
separate again as E± is decreased away from zero. Below, 
we show that Eq. 1 is only approximately true. 

In order to compare the theory of spin-splitting to the 
experimental results, we have performed self-consistent 
subband calculations that have no adjustable parameters 
(further details are in Refs. Q and Q). These calcula- 
tions produce, for a given 2D hole density and E±, both 
the B = spin-subband densities and the SdH oscilla- 
tions. It is important to note that both of these results 
are produced by the same Hamiltonian, so they are di- 
rectly comparable. These calculations yield a surprising 
result: FFTs of the calculated SdH oscillations show that 
the frequencies present, when converted to spin-subband 
densities using Eq. 1, do not agree with the calculated 
B = spin-subband densities [4|. 

This is highlighted in Fig. i which shows the calcu- 
lated B — spin-subband densities (right axis), and the 
peak positions of FFTs (fsdn) of both the calculated 
and measured SdH oscillations (left axis), as a function 
of E± |2q] . The two y-axes of the figure are scaled by 
Eq. 1, so the sets of data can be directly compared. The 
calculated fsdH consistently underestimate the B = 
spin-splitting. The agreement between the theoretical 
and measured fsdH, however, is very good. Especially 
noteworthy are the results for E± ~ 0. As the B = cal- 
culations in Fig. |i| indicate, even at E± = 0, there should 
be finite spin-splitting because of the inversion asymme- 
try of the GaAs (zincblende) crystal structure. However, 
both the measured and calculated FFTs show only one 
peak for -1 < E ± < 1 kV/cm ||]. 

We believe that the inaccuracy of Eq. 1 is due to a 
breakdown of Onsager's argument |27|] , which is based 
on Bohr-Sommerfeld quantization of the semiclassical 
motion of the Bloch carriers. The presence of spin- 
orbit interaction makes the system inherently quantum- 
mechanical, so a semiclassical picture fails. This finding 
is quite general in that one can expect deviations from 
a semiclassical picture for all systems with strong spin- 
orbit interaction. However, the full quantum-mechanical 
calculations predict that Eq. 1 can be quite accurate for 
some 2D systems grown on high-symmetry crystal direc- 
tions, while being inaccurate for other crystal directions 
or other systems B . 
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IV. METALLIC BEHAVIOR AND 
SPIN-SPLITTING 



For many years, it was widely accepted that there can 
be no metallic phase in a 2D carrier system [^8|,^9| . How- 
ever, recent experiments on several different high quality 
2D systems have provided us with reason to re-visit this 
belief, as they showed that p has a metallic-like tem- 
perature dependence: at very low temperatures p de- 
creases with decreasing T |q-|l3[]. Various mechanisms 
pi| [$H have been proposed to explain the metallic be- 
havior, but no clear model has emerged which quantita- 
tively describes the sizeable body of experimental data. 

By tuning the spin-splitting at constant density, we 
isolate the effect of spin-splitting on the T-dependence 
of p (l^Jl^]. Figure ||C shows p as a function of T for 
each of the measured E±. The traces are lined up with 
the corresponding FFTs in Fig. ||B. From this data it is 
evident that the magnitude of the change in p from 30 
mK to ^0.7 K is correlated with the spin-splitting. The 
30 mK, B — resistivity (po) for each trace is listed on 
the y-axis. 

In order to characterize the T-dependence data in a 
simple way, we calculate Ap T /po, the fractional change 
in resistivity from 30 mK to ~0.7 K. This is plotted in 
Fig. H vs. the spin-subband population difference cal- 
culated at B = (Ap s ). We have repeated this type 
of experiment on lower density samples, down to a den- 
sity of 2.5 x 10 10 cm~ 2 . For densities from 7.0 x 10 10 to 
3.3 x 10 11 cm~ 2 , for which the p vs. T traces have quali- 
tatively the same shape as in Fig. |||C, we have calculated 
the B = spin-subband densities. As shown in Fig. ||, 
data for these densities exhibit the same trend: Ap T / po 
is larger at larger Ap s , and Ap T / po is more sensitive to 
spin-splitting at the lower measured densities. 

The T-dependence of p for a density of 2.5 x 10 10 
cm~ 2 has a qualitatively different shape from the higher- 
density data: it exhibits a local maximum. Increasing 
Ap s moves this local maximum to lower T. We note that 
previous experiments on the density-dependence of the 
metallic behavior show that at high densities p monoton- 
ically increases with T, and that as the density is reduced 
and the transition to insulating behavior is approached, 
a local maximum appears in p vs. T [H||-0. Qualita- 
tively, in the density range we have measured, the effect 
of increasing spin-splitting on shape of p vs. T is similar 
to the effect of reducing density. 

For the highest density, 3.3 x 10 11 cm -2 , we have di- 
rect experimental support for the calculated Ap s . At the 
lower densities the sample quality is typically worse and 
the spin-splitting is too small so that two frequencies are 
not resolved in the SdH oscillations. We use the calcula- 



tions described in Section III to determine the expected 
Ap s used in Fig. || from the E±. Even if a large er- 
ror in Ap s is allowed for, the conclusions of the previous 



paragraphs are still valid. 

Tuning E± tunes the spin-splitting, but can also affect 
the mobility and cause changes in pq. However, a careful 
examination of all our data reveals that the changes in po 
are not causing the changes in Ap T / p$. The variation of 
Ap T / po due to E± does not correlate with the changes 
in po- One example of this is in Fig. ^C, where there are 
three traces at different E±_ that have different Ap T /po, 
but the same po (58.3 fi/sq.). 

All of the data we have presented so far have been from 
the low-mobility [011] arm of the T-shaped Hall bar. The 
SdH oscillations in the high-mobility [233] direction data 
are very similar to those in the [011] traces, and the FFTs 
show that the frequencies present, as expected, are the 
same. The T-dependence data, while qualitatively sim- 
ilar, are different along [233]. In a given measurement, 
Ap T /po is always smaller for the [233] direction than for 
the [011] direction. An example of this behavior can be 
seen in Fig. ||B. 

In summary, we have found that in the density regime 
where p increases monotonically with increasing T, the 
magnitude of the change in p, Ap T /po, is correlated with 
the spin-splitting of the 2D hole system. As the density 
is reduced, Ap T /p becomes larger and more sensitive to 
increased spin-splitting. In the density regime where p 
has a local maximum, increasing the spin-splitting moves 
the maximum to lower T. We also find that the direction 
of the current in the sample plays a surprising role: the 
higher-mobility direction shows a smaller Ap T / p$. 



V. IN-PLANE MAGNETIC FIELD 

We also employed an in-plane magnetic field, using the 
tilting stage of the 3 He refrigerator, to probe the Zeeman 
splitting in GaAs 2D holes (34]. Similar magnetoresis- 
tance (MR) measurements have been recently reported 
for 2D systems that exhibit a low-T metallic behavior 
|Io| , ^4| -[T8i1 . In our measurements we have discovered a 
remarkable anisotropy in the effect of B on the spin- 
subbands, pointing out that the Zeeman splitting in this 
2D system is very anisotropic. 

The measurements were done on a sample similar to 
those described above, but the densities were lower still. 
As with previous samples, we used both front and back 
gates to control the density. The longitudinal and Hall 
MRs were first measured with a perpendicular £?, and 
then the sample was tilted 90° and p as a function of 
in-plane B was measured. The measurements were made 
once with the sample mounted with the [233] direction 
parallel to the tilt axis. Then the sample was warmed 
up, re-mounted with the [011] axis parallel to the tilt 
axis, and the measurements were repeated. Note that 
this sample also had an T-shaped Hall bar geometry as 
shown in Fig. ^B, so for each orientation of the sam- 
ple relative to the in-plane B, p was measured for both 
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current directions simultaneously. 

In order to make the data easier to track, we scale the 
traces to the B = 0, base-T value of p. Figure ^ shows 
such data for various densities, organized by the relative 
orientations of B and the crystal axes. All traces have 
an overall positive MR and, in addition, show a broad 
feature: there is an inflection point followed by a reduc- 
tion in slope, followed by another inflection point beyond 
which the traces curve upwards again. To highlight this 
behavior, the arrows in Fig. ^ are placed between the 
two inflection points, at a B we will refer to as B* . Sur- 
prisingly, for each density, B* for the B || [233] traces 
is about 4 T smaller than for the B || [Oil] traces, re- 
gardless of the I direction. Also, B* becomes smaller as 
the density is reduced. Figure || reveals that the relative 
orientations of B and the crystal axes play an important 
role in the MR features. 

The existence of the MR features around B* is in- 
triguing. Similar, though sharper, features have been 
observed in in-plane B measurements in systems with 
multiple confinement subbands when a subband is de- 
populated ]35|] . Similarly, the features in our data may 
be related to a spin-subband depopulation and the result- 
ing changes in subband mobility and inter-subband scat- 
tering as the in-plane B is increased. In support of this 
hypothesis, numerical calculations similar to the ones de- 
scribed above show that the spin-subband depopulation 
happens at much lower B for B || [233] than for B || [Oil] 
p4j . This reflects the highly anisotropic spin-subband 
structure of the 2D hole systems in GaAs (311)A QWs. 
Our hypothesis is further supported by the observation 
(Fig. |6|) that while B* clearly depends on the orientation 
of B with respect to the crystal axes, it is rather insensi- 
tive to the current direction: spin-subband depopulation 
should not depend on the direction of current in the sam- 
ple. 

At higher in-plane B, beyond the MR features around 
B* , the data in Fig. || are qualitatively similar. The 
traces for B _L / have greater slope than the correspond- 
ing traces with B \\ I, regardless of crystal axes. In this 
regime the magnetic confinement can become compara- 
ble to the electric confinement, and the effects due to the 
finite-thickness of the 2DHS may be dominant. Indeed, 
Ref. |2^] predicts that MR with in-plane B should be sig- 
nificantly larger for B _L / than for B || J, in agreement 
with our highest B data. 

Figure shows the T-dependence of the MR at a den- 
sity of 3.9 x 10 10 cm~ 2 , for the four measured relative 
orientations of B, I, and the crystal axes. For each 
panel, the traces exhibit a nearly T-independent mag- 
netic field Bt which occurs near the trace's first inflec- 
tion point. This is consistent with the data of Ref. Jl8| . 
For B < Bt, the data show metallic behavior, and for 
B > Bt, insulating behavior. Bt is different in each 
panel and, similar to B* , depends much more strongly 
on the orientation of the crystal axes relative to B than 



on the orientation of / relative to B. Our experiments 
indicate that B* and Bt depend very similarly on the 
parameters of our systems {p, E±, direction of B). Our 
observation, which is in agreement with the in-plane MR 
data of Ref. |l(| , strongly suggests that the metallic be- 
havior is linked to the presence of two populated spin- 
subbands [[2|Jl3|jl|,|2§ . 



VI. SUMMARY 

We have demonstrated tunable B — spin-splitting at 
constant density in a GaAs 2D hole system. In the pro- 
cess, we have discovered that systems with a significant 
spin-orbit interaction show a more complicated relation- 
ship between the B = spin-subband populations and 
the frequencies present in the Shubnikov-de Haas oscil- 
lations than had previously been expected. Using the 
tunability of the spin-splitting to investigate its effect on 
the metallic behavior observed in this 2D system, we find 
that changing the magnitude of the spin-subband popu- 
lation difference changes the T-dependence of p. 

Through the use of an in-plane magnetic field, we have 
measured a surprising anisotropy of the subband struc- 
ture of (311)A GaAs 2D holes. When a magnetic field 
is applied in the plane of the 2D system parallel to the 
[233] direction, the upper spin-subband is depopulated 
at a significantly lower field than if the field is applied 
parallel to the [Oil] direction. Furthermore, we observe 
that the B = metallic behavior turns into insulating 
near B at which the upper spin-subband depopulates. 

Finally, we note that Das Sarma and Hwang have re- 
cently reported calculations aiming to explain the T- 
dependence of the resistivity J3j| and the in-plane MR 
fl22| of 2D systems that exhibit metallic behavior at finite 
T. Their calculations, which include only charged impu- 
rity scattering and the orbital motion, qualitatively re- 
produce some of the experimental data. We wish to point 
out that our results reveal the importance of the spin de- 
gree of freedom, and suggest that for an understanding of 
the experimental data it is important to also consider a 
scattering mechanism involving the spin-subbands, per- 
haps intersubband scattering [|l9|j2(J . Also important for 
(311)A GaAs 2D holes is the inclusion of interface rough- 
ness scattering: both the T-dependence of p at B = 
(Fig. IJd), as well as the in-plane MR data (Fig. ||), 
depend on the direction of the current in the crystal. 

This work was supported by the ARO and the NSF. 
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FIG. 1. A: Schematic cross section of the sample, not to scale. B: Diagram of the L-shaped Hall bar used for measuring the 
resistivities along the [Oil] and [233] directions. C: Schematic demonstrating how front and back gates can be used to tune the 
symmetry of the quantum well, and therefore the spin-splitting, without changing the density. 
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FIG. 2. A: Resistivity p data for magnetic field B perpendicular to the plane of the 2D hole system with a density p = 6.3 x 10 
cm~ 2 and current I | [233] at T — 30 mK. The data exhibits fractional quantum Hall effect at low filling factors (y), 
demonstrating the high quality of the sample. The inset shows low-B data for I || [011] (solid trace) and / || [233] (dashed 
trace). B: B = temperature-dependence at p = 3.3 x 10 10 cm -2 , highlighting the difference between [011] (solid) and [233] 
(dashed) directions. The main figure shows the fractional change in p as T is increased, while the inset shows the raw data. 
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FIG. 3. A: Magnetoresistance traces, all at a density of 3.3 x 10 11 cm~ 2 , but at different values of E±. The data shown 
are from the low-mobility [Oil] direction. B: Fourier transforms of the Shubnikov-de Haas oscillations, showing that the 
spin-splitting is being tuned through a minimum. C: Temperature dependence of p for the [Oil] direction. The traces are 
shifted vertically for clarity, with the value of p at 30 mK listed along the j/-axis for each trace. Each p vs. T trace is aligned 
with its corresponding Fourier transform. B and C together show that the magnitudes of the spin-splitting and the temperature 
dependence are related. 
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FIG. 4. Calculated B = spin-subband densities, and the peak positions of FFTs of both the calculated and measured SdH 
oscillations. The left and right axes are related according to Eq. 1 so the data can be directly compared. 
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FIG. 5. Fractional change in resistivity, measured along [011], from 30 mK to ~0.7 K plotted vs. the calculated spin-subband 
population difference. It is evident that the magnitude of the T-dependence becomes more sensitive to spin-splitting at lower 
densities. 
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FIG. 6. Fractional change in resistivity due to an in-plane B, showing that the relative orientations of B and the crystal axes 
play an important role in determining the position of the magnetoresistance features. The vertical arrows mark B* as defined 
in the text. 
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FIG. 7. Magnetoresistance data at various temperatures, for density p = 3.9 x 10 10 cm 2 , for the four relative orientations 
of B, I, and crystal axes. The fields Bt at which the resistivity is nearly T-independent are indicated by vertical marks. 
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